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PERCOLATION GAMES, PROBABILISTIC CELLULAR 
AUTOMATA, AND THE HARD-CORE MODEL 


ALEXANDER E. HOLROYD, IRENE MARCOVICI, AND JAMES B. MARTIN 


Abstract. Let each site of the square lattice Z 2 be independently declared 
closed with probability p , and otherwise open. Consider the following game: a 
token starts at the origin, and the two players take turns to move it from its 
current site x to an open site in {x + (0,1), x + (1, 0)}; if both these sites are 
closed, then the player to move loses the game. Is there positive probability 
that the game is drawn with best play i.e. that neither player can force a 
win? This is equivalent to the question of ergodicity of a certain elementary 
one-dimensional probabilistic cellular automaton (PCA), which has been studied 
in the contexts of enumeration of directed animals, the golden-mean subshift, 
and the hard-core model. Ergodicity of the PCA has been noted as an open 
problem by several authors. We prove that the PCA is ergodic for all 0 < p < 1, 
and correspondingly that the game on Z 2 has no draws. We establish similar 
results for a certain misere variant of the game and a PCA associated to it. 

On the other hand, we prove that the analogous game does exhibit draws for 
sufficiently small p on various directed graphs in higher dimensions, including 
an oriented version of the even sublattice of Z d in all d > 3. This is proved via a 
dimension reduction to a hard-core lattice gas in dimension d— 1. We show that 
draws occur whenever the corresponding hard-core model has multiple Gibbs 
distributions. We conjecture that draws occur also on the standard oriented 
lattice Z d for d > 3, but here our method encounters a fundamental obstacle. 


1. Introduction 

We introduce and study percolation games on various graphs. For the lattice 
Z 2 , we show that the probability of a draw is 0; this is equivalent to showing the 
ergodicity of a certain probabilistic cellular automaton. In higher dimensions, 
we prove that draws can occur, by developing a connection to the question of 
multiplicity of Gibbs distributions for the hard-core model. 
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Figure 1. Outcomes of the percolation game on the finite region 
(x G Z 2 : x\ + X 2 < 200}, declaring a draw if the token reaches the 
diagonal x± + X 2 = 200, and with p — 0.1 (left) and p = 0.2 (right). 
Colours indicate the outcome when the game is started from that 
site: first player win (blue); first player loss (green); draw (red). 
Closed sites are black. 


1.1. Two dimensional games and ergodicity. Let each site of Z 2 be closed 
with probability p and open with probability 1 — p, independently for different 
sites. Consider the following two-player game. A token starts at the origin. The 
players move alternately; if the token is currently at x, a move consists of moving 
it to x + (0,1) or to x + (1, 0). The token is only allowed to move to an open site. 
In accordance with the “normal play rule” of combinatorial game theory, if both 
the sites x + (0,1) and x + (1, 0) are closed, then the player whose turn it is to 
move loses the game. The entire configuration of open and closed sites is known 
to both players at all times. We call this the percolation game on Z 2 . 

If p > 1 — p c , where p c is the critical probability for directed site percolation, 
then, with probability 1, only finitely many sites can be reached from the origin 
along directed paths of open sites, and so the game must end in finite time. In 
particular, one or other player must have a winning strategy. (A strategy for one 
or other player is a map that assigns a legal move, where one exists, to each vertex; 
a winning strategy is one that results in a win for that player, whatever strategy 
the other player uses.) Suppose on the other hand that p < 1 — p c ] is there now a 
positive probability that neither player has a winning strategy? In that case we 
say that the game is a draw, with the interpretation that it continues for ever 
with best play. (When p — 0 the game is clearly always a draw). 
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Figure 2. The probabilistic cellular automaton (PCA) A p . 

See Figure [l] for simulations on a finite triangular region, with draws imposed as a 
boundary condition. As the size of this region tends to oo, the probability of a 
draw starting from the origin converges to the probability of a draw on Z 2 ; the 
question is whether this limiting probability is positive for any p. 

Related questions are considered in [ HMJ and [ BHMW15 J. in which the underlying 
graph is respectively a Galton-Watson tree, and a random subset of the lattice 
with undirected moves. 

In our case of a random subset of Z 2 with directed moves, the outcome (win, 
loss, draw) of the game started from each site can be interpreted in terms of 
the evolution of a certain one-dimensional discrete-time probabilistic cellular 
automaton (PCA); the state of the PCA at a given time relates to the outcomes 
associated to the sites on a given Northwest-Southeast diagonal of Z 2 . 

The PCA has alphabet {0,1} and universe Z, so that a configuration at a given 
time is an element of {0, 1} Z . (The three game outcomes will correspond to the 
two states of the PCA via a coupling of two copies of the PCA). The evolution of 
the PCA is as follows. Given a configuration r/ t at some time t, the configuration 
r/t+i at time t + 1 is obtained by updating each site n G Z simultaneously and 
independently, according to the following rule. 

• If r} t (n — 1) = rjt(n) = 0, then r] t+ \ (n) is set to 0 with probability p and 1 
with probability 1 — p. 

• Otherwise (i.e. if at least one of p t [n — 1) and r/ t (n) is 1), r] t+ i(n) is set to 
0 with probability 1. 

We denote this PCA A p , and call it the hard-core PCA. Its evolution rule at 
each site is illustrated in Figure [2] (The time coordinate t increases from top to 
bottom, and the spatial coordinate n increases from left to right). 

The hard-core PCA A p has already been studied from a number of different 
perspectives. It is closely related to the enumeration of directed lattice animals, 
which are classical objects in combinatorics. The link was originally made by Dhar 
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|Dha83] . and subsequent work includes [ BM98L 1.B.M07 ~ see also Section 4.2 of 
the survey of Mairesse and Marcovici [MM14aj for a short introduction, ft also has 
strong connections to the hard-core lattice gas model in statistical physics (which 
also has applications for example to the modeling of communications networks) - 
see Section [3] of this paper. The case p — 1/2 in particular relates to the measure 
of maximal entropy of the golden mean subshift in dynamical systems - see [ Elo96 ] 
and also [ Marl3 , Chapter 8]. 


Formally, we take A p to be the operator on the set of distributions on {0,f} z 
representing the action of the PCA; if p is the distribution of a configuration in 
{0,1} Z , then A p p is the distribution of the configuration obtained by performing 
one update step of the PCA. A stationary distribution of a PCA F is a 
distribution p such that Fp = p. (More generally, p is fc-periodic if F k p = p , 
and periodic if it is fc-periodic for some k > 1.) A PCA is said to be ergodic if it 
has a unique stationary distribution and if from any initial distribution, the iterates 
of the PCA converge to that stationary distribution (in the sense of convergence 
in distribution with respect to the product topology). 


It is straightforward to show that if p is sufficiently large, then the hard-core 
PCA A p is ergodic. The question of whether A p is ergodic for all p G (0, f) has 
been mentioned as an open problem by several authors, for example in |TVS + 90l| . 
[ LBM07 ]. [ MM14a ], A link between the hard-core PCA A p and the percolation 
game was already mentioned by [LBM07] , and it is relatively easy to show that the 
percolation game has positive probability of a draw if and only if A p is non-ergodic. 


We also consider a certain “misere” variant of the game which we call the target 
game. In the percolation game described above, moves to closed sites were 
forbidden; this is of course equivalent to allowing all moves, but with the proviso 
that a player who moves to a closed site loses the game. In the target game, the 
game is instead won by the first player who moves to a closed site. In the same 
way that the percolation game corresponds to the hard-core PCA A p , we show 
that the target game corresponds to another PCA, which we denote by B p , whose 
update rule is shown in Figure [3j 


PCA that are defined on Z and whose alphabet and neighbourhood are both of size 
2 are sometimes called elementary PCA. A variety of tools have been developed to 
study their ergodicity. Under the additional assumption of left-right symmetry of 
the update rule, these PCA are defined by only three parameters: the probabilities 
to update a cell to state 1 if its neighbourhood is in state 00, ff, or Of (which 
is the same as for fO). Existing methods can be used to handle more than 90% 
of the volume of the cube [0, f] 3 defined by this parameter space, but when p is 
small, the PCA A p and B p belong to an open domain of the cube where none of 
the previously known criteria hold |TVS + 90l Chapter 7]. 
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1 with probability p 


Figure 3. The PCA B p . 


We now state our first main result. 

Theorem 1. For any p e (0,1), the PCA A p and B p are both ergodic, and the 
probability of a draw is zero for both the percolation game and the target game on 


We prove ergodicity by considering the envelope PCAs corresponding to A p and 
B p , which are PCAs with an expanded alphabet {0, ?, 1}. The envelope PCA 
corresponds to the status of the game started from each site (with the symbols 0, 
? and 1 corresponding to wins, draws and losses respectively). An evolution of 
the envelope PCA can be used to encode a coupling of two copies of the original 
PCA, with a ? symbol denoting sites where the two copies disagree. We introduce 
a new method involving a positive weight assigned to each ? symbol (whose value 
depends on the state of nearby sites). The correct choice of weights is delicate 
and non-obvious. We show that if the process is translation-invariant, then the 
average weight per site strictly decreases under the evolution of the envelope PCA, 
unless it is 0. It follows that any translation-invariant stationary distribution for 
the envelope PCA has no ? symbols, with probability 1, and from this we will be 
able to deduce that the game has no draws with probability 1, so that the original 
PCA is ergodic. Although the proof of ergodicity could be phrased so as not to 
refer to games, the notion is useful as a semantic tool and a guide to intuition. 
Our proof of Theorem fl] follows almost identical steps for the PCA B p (and the 
target game) as for A p (and the percolation game). However, in some respects the 
two cases are very different. In particular, for the percolation game, combining 
Theorem [T| with known methods permits an explicit description of the distribution 
of game outcomes along a diagonal, as a Markov chain. Consequently, we show 
that the probability that the first player wins the percolation game or the origin 
is closed is 


( 1 . 1 ) 
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P 

Figure 4. The probability that the first player wins the percolation 
game, conditional on the origin being open, as a function of the 
density p of closed sites. 


See Figure [4] for the conditional first player win probability given that the origin 
is open. This conditional probability is greater than 1/2 if and only if p £ (0,1/3), 
and is maximized at p = (2 — a/3)/3 = 0.0893.... These methods seemingly do 
not work for the target game, and we do not know an explicit expression for 
the win probability. Furthermore, for an alternative misere version of the game 
(see Question 4.5 in the final section of the paper) there is apparently no similar 
connection to a PCA with alphabet {0,1}, and we do not have a proof that the 
probability of a draw is 0. This is somewhat reminiscent of the situation for sums 
of combinatorial games [Con nu, where the well-developed theory of “normal play” 
games extends only in very limited cases to their misere cousins. 


1.2. Ffigher dimensions and the hard-core model. We also consider exten¬ 
sions of the percolation game, described above for Z 2 , to more general directed 
graphs, and in particular to lattices in higher dimensions. Theorem [l] tells us that 
in two dimensions, the probability of a draw is 0 for all positive p, but we find a 
very different picture in three and higher dimensions. 

Let G = (V, E) be a locally finite directed graph. For x G V, let Out(x) and 
In(a;) be the sets of out-neighbours and in-neighbours of x respectively. For the 
percolation game on G, let each vertex x be closed with probability p and open 
with probability 1 — p, independently for different vertices. A token starts at some 
vertex, and the two players move alternatively; if the token is currently at x, a 
move consists of moving it to any vertex in Out(x). The token is only allowed 
to move to open sites; if all the vertices in Out(x) are closed, then the player to 
move loses the game. 
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For graphs G with an appropriate structure, we develop a connection to the 
hard-core model on a related undirected graph in one fewer dimensions, to obtain 
a criterion under which the game is drawn with positive probability. 

For an undirected graph with vertex set W, and any A > 0, a Gibbs distribution 
for the hard-core model on W with activity A is a probability distribution on 
configurations r) G {0, 1} W such that 


(1.2) P^(f) — 1 (v( w ) : w / v)^j 


A 


if r](w) = 0 for all 
neighbours w of u; 


1 + A’ 

0, otherwise. 


Any such Gibbs distribution is concentrated on configurations p that correspond 
to independent sets, in the sense that no two neighbouring vertices v and w have 
i](v) = r]{w ) = 1. If W is finite, then there is a unique Gibbs distribution, which is 
the probability distribution that puts weight proportional to on eac h 

configuration r] that is supported on an independent set. However, for infinite 
graphs, there may be multiple Gibbs distributions. A well-known example is 
the lattice Z d with nearest-neighbour edges. For d — 1, there is a unique Gibbs 
distribution for all activities A. However, for d > 2, there exist multiple Gibbs 
distributions when A is sufficiently large | Dob65 j. 

Returning to the percolation game on a directed graph G, we now give the key 
assumptions on G that are required for our dimension reduction method. Suppose 
there is a partition (S'*. : k G Z) of the vertex set V of G, and an integer m > 2, 
such that the following conditions hold. 


(Al) For all x G S k , we have Out(a;) C S k +i U • ■ ■ U S k + m - i. 

(A2) There is a graph automorphism 0 of G that maps S k to Sk+ m for every k , 
and such that Out(x) = ln(0(a;)) for all x. 


Then let Dk be the graph with vertex set S'*, U ■ ■ • U Sk+ m - 1 , with an undirected 
edge ( x,y ) whenever (x, y) is a (directed) edge of V. (Below for convenience we 
will also use D k to denote the vertex set S k U • • • U Sk+ m - 1 -) It is straightforward 
to show that under conditions (Al) and (A2), the graphs D k are isomorphic to 
each other for all k G Z (see Lemma 3.1); write D for a generic graph which is 
isomorphic to any of the D k . Then we have the following criterion for positive 
probability of draws. 


Theorem 2. Suppose that the directed graph G satisfies (Al) and (A2). If there 
exist multiple Gibbs distributions for the hard-core model on D with activity X, 
then the percolation game on G with p = 1/(1 -F A) has positive probability of a 
draw from some vertex. 







HOLROYD, MARCOVICI, AND MARTIN 


The simplest case in which to understand the conditions (Al) and (A2) is when 
G is the directed latti ce Z 2 , with Out(x) = {x + (1, 0), x + (0,1)} (the setting of 
Theorem [l] in Section 1.1). Then we may take the partition of Z 2 into Northeast- 


Southwest diagonals given by Sk '■= {(or,£ 2 ) ■ x± + X 2 = k}, along with the 
bijection 0(x) = x + (1,1), and m = 2. The graph D then consists of the vertices 
of two successive diagonals, and is thus isomorphic to the line Z. (In the context 
of PCA, D is sometimes called the doubling graph.) 


Since, as noted above, there is a unique Gibbs distribution for the hard-core model 
on Z for all A > 0, Theorem [2] does not imply the existence of draws for any 
p € (0,1). (Indeed, that would contradict Theorem [I]). 


In higher dimensions the picture is different. We will give several examples of 
relevant graphs in Section T2 and Theorem [3] below. For the current discussion, 
consider the case where G has vertex set Zf ven := (x G Z d : Yh x i is even}, with 
directed edges given by Out(x) := {x ± e* + : 1 < i < d — 1} (where e* is 
the zth standard basis vector in Z d ). So Out(x) has size 2{d — 1); any move 
of the game increases the cith coordinate by 1 and also changes exactly one of 
the other coordinates by 1 in either direction. In two dimensions, this game is 
isomorphic to the original game on Z 2 . For general d, conditions (Al) and (A2) 
hold with m — 2 if we set Sk = {x G Zf ven : Xd = k} and </>(x) = x + 2e^. One 
then finds that D is isomorphic to the standard (d — l)-dimensional cubic lattice 
Z d_1 with nearest-neighbour edges. As mentioned above, there are multiple Gibbs 
distributions for the hard-core model on Z d_1 whenever d > 3 and A is large 
enough; then Theorem [2] tells us that the percolation game on G has positive 
probability of a draw when p is sufficiently small. We do not know whether the 
draw probability is monotone in p. nor even whether it is supported on a single 
interval (giving a single critical point). 


To prove Theorem [2j we consider a recursion, analogous to the earlier PCA, 
expressing game outcomes starting from vertices in Sk in terms of outcomes 
starting in Sk+i U • • • U Sk + m ■ Via the graph isomorphism from Dk to D, the 
iteration of this recursion can be reinterpreted as a version of Glauber dynamics for 
the hard-core model on D. If the hard-core model has multiple Gibbs distributions, 
then they correspond to multiple stationary distributions for the recursion on G, 
and from this we will deduce that draws occur. 


The stationary distributions for the recursion on G that arise from the above 
correspondence have a certain special property, which in the case m = 2 can be 
viewed as a version of time-reversibility. I 11 cases where the hard-core model has 
a unique Gibbs distribution (such as when G — Z 2 and D — Z), the argument 
implies that the recursion has only one stationary distribution with this reversibility 
property, but says nothing about the possibility of other stationary distributions. 
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It is for this reason that the implication in Theorem [2] is in only one direction, and 
that the proof of Theorem [l] requires a different argument. 

Unfortunately, the following very natural example is not amenable to our methods. 
Let G be the standard cubic lattice 7L d with edge orientations given by Out(x) = 
{x + e, : 1 < i < d}. Theorem [2] does not apply for d > 3, because there is no 
choice of m and the automorphism 0 such that (A2) holds. We conjecture that, 
nonetheless, the percolation game has positive probability of a draw whenever p is 
sufficiently small. 


1.3. Further background. The celebrated positive rates conjecture is the asser¬ 
tion that in one dimension, any finite-state finite-range PCA is ergodic, provided 
the transition probability to any state given any neighbourhood states is positive 
(the latter is the “positive rates” condition). This contrasts with two and higher 
dimensions, where for example the low temperature Ising model is well-known to be 
non-ergodic. Despite persuasive heuristic arguments in favour of the positive rates 
conjecture, Gacs jGacOlj has given an extremely complicated one-dimensional 
PCA refuting it. (See also j GraOlj .) However, it is still natural to hypothesize 
that all “sufficiently simple” one-dimensional PCA with positive rates are ergodic. 


The PCA A p and B p are very simple, but do not satisfy the positive rates condition. 
(E.g. in Ap , the word 11 deterministically yields 0.). However, similar but weaker 
conditions do hold; for example, any finite word in {0, l} n has positive probability 
of yielding any word in {0, l} n ~ 2 after two steps of the evolution. In light of this 
and the above remarks, it would have been very surprising if these PCA were not 
ergodic. Nonetheless, proving ergodicity is often very difficult, even in cases where 
it appears clear from heuristics or simulations. 

Another case in point is the notorious noisy majority model on 7L d . Here, a 
configuration is an element of {0, 1} Z . The update rule is that with probability 
1 — p, a site adopts the more popular value in {0,1} among itself and its 2d 
neighbours; with probability p it adopts the other value. In dimensions d > 2 it 
is expected that this PCA should behave similarly to the Ising model: it should 
be ergodic for p sufficiently close to 1/2, and non-ergodic for p sufficiently small, 
with a unique critical point separating the two regimes. However, proving any 
of this appears very challenging. See e.g. [BBJW10 1 IGraOl] and the references 
therein for more information. One key difficulty with the noisy majority model is 
the lack of reversibility of the dynamics (in contrast to the Glauber dynamics for 
the Ising model, for example). This can be compared to the difficulty of obtaining 
a result like Theorem [2] in the absence of conditions such as (Al) and (A2); see 
the discussion above at the end of Subsection O 


In a different direction, a variant of the notion of envelope cellular automata has 
recently been combined with percolation ideas in jGHISj . to prove the surprising 
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fact that certain deterministic one-dimensional cellular automata exhibit order 
from typical finitely supported initial conditions, but disorder from exceptional 
initial conditions. 


1.4. Organization of the paper. In Section [2] we explain the link between the 
PCA A p and B p and the percolation game and target game respectively in Z 2 . We 
also establish several basic results concerning monotonicity and ergodicity. The 
local weighting on configurations is introduced in Subsection 2.3, and the proof of 


ergodicity is then given in Subsection 2.4 


The relation to the hard-core model is then developed in Section [3j We start by 
considering the case of Z 2 where the ideas are simplest, and in particular we will 
derive the formula (1.1) for the winning probability. The general case is then 
treated in Subsection T2 where Theorem [2] is proved. In Subsection |3.4| and 
Theorem [3j we give a variety of examples of the application of Theorem [2] to 
graphs with vertex set 7L d for d > 3, for which the role of the doubling graph D is 
played by various lattice structures. We also give an extension of Theorem [2] in 
Proposition 3T in Subsection 3.5 using a variant form of the correspondence to 
the hard-core model. 

We conclude in Section [4] with some open problems. 


2. Percolation games and probabilistic cellular automata 

2.1. The PCA for the percolation game. Consider the percolation game on 
Z 2 as defined in the introduction. 

Suppose x is an open site of Z 2 . Let rj(x) be W, L or D according to whether the 
game started with the token at x is win for the first player, a loss for the first 
player, or a draw, respectively. (Recall that we assume optimal play, with the 
players able to see entire configuration of open and closed sites when deciding on 
their strategies). If x is a closed site, it is convenient to set r}(x) = W. (We can 
imagine that a player is allowed to move the token to x, but with the effect that 
the game is then declared an immediate win for the opponent). 

Recall that Out (a;) = {x + ei, x + 62 } is the set of sites to which the token can 
move from x. By considering the first move, we have the following recursion for 
the status of the sites: 

x closed =>■ rj(x) = W; 

P 

=> Tj(x) = < W 

P 


( 2 . 1 ) 


x open 


if rj(y) = W for all y € Out(x) 
if rj(y) = L for some y e Out(x) 
otherwise. 
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* 1 


0 0 1 * 



0 w. prob. p 0 

1 w. prob. 1 — p 


? 0 

0 ? 

? ? 



0 w. prob. p 

? w. prob. 1 — p 


Figure 5. The PCA F p (* denotes an arbitrary symbol). 


For k G Z, let Sk be the s et {x = (x i, x 2 ) G Z 2 : Xi + x 2 = k}, a NW-SE diagonal 
of Z 2 . The recursion (2.1) gives us the values (r)(x) : x G Sk) in terms of the values 
: x G Sk+i) together with the information about which sites in Sk are closed. 

It is important to note that it is not a priori clear whether the recursion (2.1) 
suffices to determine rj uniquely from the configuration of open and closed sites. 
Indeed, in the trivial case p — 0 when all sites are open, we have p{x) = D for all 
x, but (2.1) has other solutions: one is to set rj'(x) equal to L or W according to 
whether x\ + x 2 is odd or even. Such considerations are in fact central to many of 
our arguments. One way to interpret our main result, Theorem [T] is as saying that 
(2.1) does have a unique solution almost surely for all 0 < p < 1 (and similarly 
for the target game). In contrast, for the higher dimensional variants considered 
later, the analogous recursions admit multiple solutions for suitable p. 

Via (2.1), we can regard the configurations on successive diagonals Sk, as k 
decreases, as successive states of a one-dimensional PCA. Let us introduce the 
following recoding: 

W = 0, L = 1, D = ?. 


(In the coupling arguments below, the symbol ? will be interpreted as marking 
a site at which the value is “unknown”. The choice to assign W = 0 and L = 1, 
rather than the other way around, say, will be important for the later connection 
with hard-core models.) The PCA evolves as follows: given the values for sites in 
Sk- i-i, each value rj(x) for x G Sk is derived independently using the values r)(x + e i) 
and p(x + e 2 ), according to the scheme given in Figure [5] (where a * represents an 
arbitrary symbol in {0, ?, 1}). 

We denote the corresponding PCA F p . Although we have defined it as a process in 
the plane, we can also regard it as a PCA on Z with a configuration in {0, ?, 1} Z 
evolving in time by setting 


(2.2) 


rit(n) = i ?((-£ - n,n)). 
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(Here we have made the arbitrary choice to offset leftward as time increases,so 


that the PCA rule gives r] t+ \{n) in terms of rj t (n) and r] t (n + 1).) As in Section 1.1 


formally we take F p to be an operator on the set of distributions on {0,?, 1}' 
representing the action of the PCA. 


In the setting of the percolation game, translation invariance of the whole process 
on Z 2 implies that the distribution of the configuration on the diagonal S'*, does 
not depend on k ; that is, the distribution of (rj((k — n,n)) : n G Z) does not- 
depend on k and is a stationary distribution of F p . In addition, this distribution 
is itself invariant under the action of translations of Z. 


We next note two useful monotonicity properties for the PCA F p . In terms of the 
game, they have natural interpretations: (i) an advantage for one player translates 
to a disadvantage for the other; and (ii) declaring draws at some positions can 
only result in more draws elsewhere. 


Lemma 2.1. Let fi and v be probability distributions on {0, ?, 1} Z . 


(i) If h A v, where < denotes stochastic domination with respect to the 
coordinatewise partial order induced by 0 < ? < 1, then F p p > F p v. (Note 
the reversal of the inequality). 

(ii) If fi <] v, where < denotes stochastic domination with respect to the 
coordinatewise partial order induced by 0 <3 ? t> 1, then F p p <j F p u. 


Proof. We can use the recursion ( 2. 1[ ) to give a coupling of a single step of 
the PCA F p started from two different configurations. Suppose we fix values 
(rj(x) : x G Sk+ 1) and (rj(x) : x G Sk+ 1), in such a way that rj(x) < rj(x) for 
all x G Sfc+i (where < is the coordinatewise order on configurations induced by 
0 < ? < 1). Now use (2.1) to obtain values rj(x) and rj(x) for x G Sk, using the 
same realization of closed and open sites in Sk in each case. It is straightforward 
to check that in that case rj(x) > rj{x) for each x G Sk- Hence the operator F p is 
decreasing in the desired sense. 


Similarly, if rj(x) < rj(x) for all x G Sk+ 1, then we obtain rj(x) < rj(x) also for each 
iG 5 (;. So in this case the operator F p is increasing as desired. □ 


If we restrict the PCA F p to configurations that do not contain the symbol ?, 
we recover precisely the hard-core PCA A p defined in the introduction. In the 
terminology of Busic et al. BMM 131 . the PCA F p is the envelope PCA of A p . A 
copy of the PCA F p can be used to represent a coupling of two or more copies of 
the PCA A p , started from different initial conditions. The symbol ? represents 
a site whose value is not known, i.e. one which may differ between the different- 
copies. 
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Specifically, consider starting copies of the hard-core PCA A p from several different 
initial conditions, represented by configurations on the diagonal Sk for some fixed 
k. As in the proof of Lemma 2.1, a natural coupling is provided by the recursion 
(2.1), using the same realization of the closed and open sites on (S r : r < k). In 


particular, let k > 0 and consider three copies r/, rj and r)' , with rj and rj started 
from arbitrary initial conditions on Sk, while ryfx) = ? for all x G Sk- (So that ry 
is maximal for the ordering < in Lemma 2.1 pi)). Then we have that rj(x) <j ry{x) 
and rj(x) < ry(x) for all x G S r with r < k. This implies that if rj(x) then 

ry[x) = ?. 

In terms of the game, we have the following interpretation: if the origin O = (0, 0) 
is an open site, and ry(0) = 0 (respectively ry{0) = 1) then the first (respectively 
second) player can force a win within at most k moves of the game. 

The ergodicity of an envelope PCA implies the ergodicity of the original PCA, 
but the converse is not true in general. In our case, however, we can use the 
monotonicity property in Lemma 2.1[i) to show that the two are equivalent. 


Proposition 2.1. The PCA F p is ergodic if and only if A p is ergodic. 

Proof. It is clear from the definitions that if F p is ergodic, then A p is also ergodic. 
Conversely, suppose that A p is ergodic. Let /a be a distribution on {0, ?, 1} Z , 
and let <5o and ch the distributions concentrated on the configurations “all Os” 
and “all Is”. Then £o < At < (h, so by Lemma 2.1[ i), for k > 0 we have either 
FfSo < Ffu < F p 5i or F p 8 0 > F p /a > F p 5±, according to whether k is even or 
= A p Sq and F p 8 1 = A p 5\, and by ergodicity of A pi the latter two 


, s 

odd. But Ff8i 


sequences converge as k —* oo to the same distribution 7r, so F p n also converges 
to 7T. Thus F p is also ergodic. □ 

Proposition 2.2. For each p G [0,1], the percolation game has probability 0 of a 
draw if and only if A p is ergodic. 

Proof. If A p is ergodic then so is F p . and so the unique invariant distribution of 
F p has no ? symbols. But we know that- the distribution of the game outcomes 
along a diagonal Sk is invariant for F p . Hence with probability 1, there are no 
sites from which the game is drawn. 

For the converse, let u be a random configuration of open and closed sites on Z 2 
chosen according to the percolation measure. Consider any site x G So. If the 
game started from x is not a draw, then (since at each turn the player to move 
has only finitely many options) one player has a strategy that guarantees a win 
in fewer than N moves, where N G N is a finite random variable that depends 
on to. Consequently, if we assign any configuration of states 0,?, 1 to Sn and 


compute the resulting states on (S n : 0 < n < N) using the recursion (2.1) and 
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* 1 

1 * 



0 w. prob. 1 — p 

1 w. prob. p 


? 0 

0 ? 

? ? 



? w. prob. 1 — p 

1 w. prob. p 


Figure 6. The PCA G p . 


the configuration c o of open and closed sites, the resulting state at x is the same 
as its state for the percolation game on Z , 2 with configuration oj. 

Let 7 be the random configuration of game outcomes on Sq arising from uj. Also, 
fix a distribution v on { 0 , 1 } Z , and let 7 n be the configuration on So that results 
from assigning a configuration with law v to S n , independent of u>, and applying 
(2.1) as described above. By the argument in the previous paragraph, if the 
probability of a draw is 0 , then converges almost surely to 7 (in the product- 
topology). Hence also the distribution of converges to that of 7 . But 7 n has 
distribution A™v, so converges as n —> 00 to the distribution of 7 , which does 


not depend on u. Hence A p is ergodic. 


□ 


2.2. The target game and the corresponding PCA. We now turn to the 

target game, in which the winner is the first player to move to a closed site. As 
before, we can introduce a PCA to describe the status of the positions. The 


recursion ( 2 . 1 ) still applies when x is open, but now if x is closed we set instead 


rj(x) = L = 1. The equivalents of A p and F p for the target game are the PCA B p 
and G p defined via Figures [3] and [ 6 ] 

We may see the PCA B p as a composition of Stavskaya’s PCA and the flip operator. 
Stavskaya’s PCA (see for example |TVS + 90] ) is given by the local rule which sets 
Vt( n ) = 0 with probability p and otherwise rj t (n ) = max{r/ t _ 1 (n), r)t-i(n — 1)}, and 
provides a non-trivial example of a one-dimensional PCA where ergodicity fails. 
However, the ergodicity of the PCA B p has not been studied previously. As in the 
case of A p , it does not satisfy any of the previously known criteria for ergodicity. 


Using exactly the same arguments as for Proposition 2.1 and Proposition 2.2\ we 
have the following results. 


Proposition 2.3. The PCA G p is ergodic if and only if B p is ergodic. 

Proposition 2.4. The target game has probability 0 of a draw if and only if B p 
is ergodic. 
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Below we will show that, for all p > 0, any stationary distribution for F p or for 
G p is concentrated on {0, 1} Z , so that the probability of a ? symbol is 0. Using 
Propositions |2. 1| and |2.2[ or |2.3| and |2.4| respectively, we will then be able to deduce 
the ergodicity of A p and B p as required for Theorem [l] 

2.3. The weight function. We are concerned with the two PCA A p and B p on 
the alphabet {0,1}, shown in Figures [2] and [3j and their envelope PCA F p and 
G p , shown in Figures [5] and [ 6 j 

We introduce a deterministic cellular automaton D defined in Figure [7j (Note 
that D = F 0 = Go). 

We also introduce the randomization operator R® on {0, ?, 1} Z that changes each 
symbol into a 0 with probability p (making no change if it was already a 0 ), 
independently for different sites, and similarly the operator R p that changes each 
symbol into a 1 with probability p, independently for different sites. Observe that 

F p = R p o D and G p = R p o D. 

We now establish a property of the deterministic operator D. For a given configu¬ 
ration in { 0 , ?, 1} Z , let us weight the occurrences of the symbol ? as follows: 

• if a ? is followed by a 01, then it receives weight 3; 

• if a ? is followed by a 0 and then by something other than a 1 , it receives 
weight 2 ; 

• otherwise, a ? receives weight 1 . 

We say that the distribution of a configuration rj — (rji : i G Z) is shift-invariant 
if rj and (rj l+ k : % G Z) have the same distribution for each fceZ, and reflection- 
invariant if rj and (r/_j :ieZ) have the same distribution. If /i is a distribution 
and w G {0, ?, l} n is a finite word, we write ji{w) := p{rj : (r/i,..., rj n ) = w} for 
the corresponding cylinder probability. For a shift-invariant distribution /i on 
{0, ?, 1} Z , we introduce the quantity: 

M?oi) + M?0) + //(?), 
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which is the expected weight per site under ji. 

Lemma 2.2. If /a is a shift-invariant and reflection-invariant distribution on 
{0, ?, 1 } Z , then 

Dfi{7 01) + Dfi( ?0) + D/r(?) < /x(?01) + /r(?0) + //(?). 

Proof. By looking at the possible pre-images of each pattern, we obtain the 
following three equalities: 

z^(?) = M??) + Mo?) + M?o), 

D/i(?o) = M??i) + /i(o?i) + M?oi), 

D/x(?01) = 0. 

Summing, and using reflection invariance to deduce /x(0?) = /x(?0), we obtain 

(2.3) D/r(?01) + D/x(?0) + £>//(?) - //(?01) - /i(?0) 

= //(??) + A*(0?) H- At(??l) H- A«(0?1) 

<M??) + Mo?) + M?i) 

= M??) + M?o) + M?i) 

= M?)- □ 


Here is an informal way to explain the above result, which provides some insight 
into the reasons for the choice of weight. Let us consider a symmetric version 
of the weight system that we have introduced: for each symbol ?, we add its 
right-weight, as introduced above, to its left-weight, which is equal to 3 if it the 
previous letter is a 0 and if there is a 1 before it (pattern 10?), to 2 if the previous 
letter is a 0 and if there is something else than a 1 before, and to 1 otherwise. 
(Since in Lemma 2.2 we consider only reflection-invariant distributions, working 
with the symmetric weight is equivalent to working with the original weight). 


Thus, the weight of the symbol ? in the pattern 1?1 is equal to 1 -)- 1 = 2, while in 
the pattern 10??1, the weight of the first ? symbol is 3 (left) + 1 (right) = 4. and 
the weight of the second one is equal to 1 + 1 = 2. 

Figure [8] shows an example of evolution of the deterministic CA D from an 
initial configuration represented at the top (with time going down the page). 
The symmetrized weights of the symbols ? appearing in the space-time diagram 
are shown in red. As illustrated in the figure, from a pattern 1? 1, the symbol ? 
disappears and the weight thus decreases, but in other cases the total weight is 
locally preserved. 
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r /'\ 


' \ / 2 \ / 2 \ / 2 \ 

0 ? ? 0 

. / \ /3\ /3\ / 
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' \ / 2 \ / 2 \ / 2 \ 

0 ? ? 0 



3 3 



1 

1 

1 


FIGURE 8 . Example of evolution of the weight of a configuration 
under the operator D. Time runs down the page, and the weight of 
each ? symbol is given below it. 


2.4. Proof of ergodicity. 

Proposition 2.5. For 0 < p < 1, the PC A F p has no stationary distribution in 
which the symbol ? appears with positive probability. 


Proof. It suffices to show that there is no shift-invariant and reflection-symmetric 
stationary distribution in which the symbol ? appears with positive probability. For 
consider iterating the PCA starting from the distribution 5? concentrated on the 
configuration with ? at all sites. By Lemma 2.1 (ii), the probability Ff5?{f) is non¬ 
increasing, and if there is any stationary distribution /1 with positive probability 
of ?, then Ff <5?(?) is bounded below by /i(?) for all n, and so does not converge to 
0. Then any limit point of the sequence of Cesaro sums of FfS-? is a stationary 
distribution that has positive probability of ?, and that is also shift-invariant and 
reflect ion-symmetric. 


For a shift-invariant distribution v on {0, ?,1} Z , we have the following three 
equalities: 


i#/(?) = (l-p)i/(?), 

R 0 p vf 0) = p( 1 - p) z/(?) + (1 - p) 2 p(?0) 

R p u(701) = p( 1 — p) 2 z/(?*l) + (1 — p) 3 p(?01). 

Here * represents an unspecified symbol to be summed over, so that z/(?*l) : = 
£ a= o ? 1 p(?al)- Adding together these three equalities, we get 

R° p p(?) + R° p u(l 0) + R° p u(l 01) = i/(?) + p(?0) + p(?01) + p( 1 - p) 2 p(?*1) 

(2.4) — p 2 i/(?) — p (2 — p) p(?0) — p (3 — 3 p + p 2 ) z/(?01). 
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Let n be a stationary distribution of F p that is shift-invariant and reflection- 
symmetric, and set v = Dp. Then, since F p = R p o D, we have R p u = p. Hence 


by Lemma 2.2 


K?01) + ^(?0) + !/(?) < R° p u(? 01) + R° p u(? 0) + R° p u(7). 
Then, since p > 0, it follows from < \2A ) that 

(1 — p) 2 u(l* 1) > pv(l) + (2 — p)u(? 0) + (3 — 3 p + p 2 )z/(?01) 
(2.5) >pz/(?) + (2-p)z/(?0). 


We now proceed to obtain bounds for the left and right sides of (2.5) in terms 
of cylinder probabilities of p. We have p(7) = R Q p v{l) = (1 — p)u(J). Also 
z/(?0) = Dp{l 0) = /r(?01) + yu(?l) — /i(l?l) = /i(?01) + //(?1), since F p p = p and 
the pattern 1?1 has no preimage by F p . We thus obtain: 

V 


pH?) + (2 - pM?o) = V-M?) + (2 - p) K?oi) + p(?i)] 


( 2 . 6 ) 

But we have 


(2.7) 


= pH?) + (2 - p) [p(?oi) + p(?i)] + p( ? ) 


1 — p 
p 2 


> pp(? l) + (2 - p) [,,(?oi) + p(?i)] + y— p(?> 


1 — p 1 


> 2/i(?l) + /i(?01) + 


(1 - p) 2 v(l* 1) = R p v(l* 1) 




1 — p 


/!(?). 


= /l(?*l) 

= /i(?01) +/i(??l) + /r(?ll) 
< /r(?01) + 2/i(?l). 


Putting together (2.5), (2.6) and (2.7) we obtain that /i(?) = 0 as required. □ 

We now turn to the PCA corresponding to the target game. 

Proposition 2.6. For 0 < p < 1, the PCA G p has no stationary distribution in 
which the symbol ? appears with positive probability. 


Proof. As in the proof of Proposition 273 we need only consider stationary dis¬ 
tributions that are shift-invariant and reflection-symmetric. For a shift-invariant 
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distribution v on {0,?, 1 } Z , we have the following equalities: 

R^(?) = (1 - pM?) 

/fX?0) = (i-p)M?0) 

Rpisfi 01) = p( 1 — p) 2 z/(?0) + (1 — p) 3 z/(?01). 

Thus, 

i?>(?) + i?>(?0) + i?>(? 01) = i/(?) + z/(?0) + z/(?01) + p(l - p) 2 z/(?0) 

— pz/(?) — p( 2 — p)z/(?0) — p(3 — 3p + p 2 )z/(?01). 


Let /i be a shift-invariant, reflection-invariant, invariant distribution of G p , and 
let v = Dp. Then, R p u = p. By Lemma 

i/(?01) + p(?0) + i/(?) < i?Ji/(?01) + i?>(?0) + i?X?). 

It follows that 


2.2 


(1 — p) 2 z/(?0) > z/(?) + (2 — p)z/(?0) + (3 — 3p + p 2 )z/(?01) 

> (2 — p)p(?0), 

which is possible only if p(? 0) = 0. We then obtain /r(?0) = 0, and it follows easily 
that /r(?) = 0 . □ 


Now we can quickly deduce our main result. 


Proof of Theorem [Ij We know that the distribution of the states (win, loss, draw) 
of the sites along a diagonal Sk in the percolation game is a stationary distribution 


for F p . Since by Proposition 2.5, F p has no stationary distribution with positive 


probability of ? for any p > 0 , the probability of a draw in the percolation game 
must be 0. Then by Proposition 2.2 the PCA A p is ergodic for each p > 0. 

An identical argument applies for the target game, and the PCA B p , using 
Propositions |2.6| and |2.4[ □ 


3. Percolation games and the hard-core model 


3.1. The two-dimensional case. In this section we develop the relationship 
between the percolation game and the hard-core model. We start in the setting of 
Z 2 where the ideas are easiest to understand, but our main application will be in 
Section T2, when we establish a more general framework, and apply it to show 
that certain higher-dimensional games have positive probability of a draw when p 
is sufficiently small. 


Consider the hard-core PCA A p . This PCA is known to belong to a family of 
one-dimensional PCA having a stationary distribution that is itself a stationary 
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Markov chain indexed by Z [ BGM69 , lTVS + 90l lMM14bj . This distribution, )i p 
say, is the law of the stationary Markov chain on Z with transition matrix 

^2-p-y / p(4-3p) 2p 2 -3p+^p(4-3p) ^ 

2(1—p) 2 2(1—p) 2 


(3.1) 


P = 


P 0,0 
Pl.O 


Po.l 

Pl,l 


-P+Vp(4—3p) 2—p—yJ p(4—3p) 

2(1—p) 2(1—p) 


/ 


on state space {0,1}. (See Section 4.2 of |MM14a| - note that p there corresponds 
to our 1 — p). In fact, the evolution of the PC A started from y p is time-reversible 
- the distribution of the two-dimensional space-time diagram obtained (via the 


correspondence at (2.2)) is invariant under reflection in the line x\ + x 2 = k for 


any k. (In addition, the distribution y p is itself reversible as a Markov chain on Z, 
which corresponds to symmetry of the two-dimensional picture under reflection in 
the line x\ — x 2 ). 

By Theorem [lj we know that y p is in fact the unique stationary distribution of 
F p , and therefore the probability that the first player wins the percolation game 


starting from the origin is, as claimed in (1.1), 


Pp{ 0 ) 


Pl -° =bi + 

P 1,0 + Po,i 2 


P 


4 — 3 p 


(The conditional probability of a win given that the origin is open is then given by 

(M°) -p)/(i-p)-) 

An illuminating way to understand the presence of this Markovian reversible 
stationary distribution is to consider the doubling graph of the PCA, corresponding 
to two consecutive times of its evolution [ Vas78l KV80 . lTVS + 9n] . This is an 
undirected bipartite graph, connecting sites between which there is an influence 
induced by the rules of the PCA. 

As in Section 2.1[ we can think of a configuration of the PCA as indexed by 
a diagonal S*. = {(xi,x 2 ) ■ Xi + x 2 = k} of Z 2 . A time-step of the PCA then 
corresponds to moving from a configuration on S/.+ 1 to a configuration on S^. 

As before, let Out(x) = {x + ei,x + e 2 } for x G S k . The elements of Out(x) lie in 
Sk- )-i, and are the sites to which the token may move from sites x\ they are the 


sites whose values appear on the right side of the recurrence (2.1) for the value 
rj(x). Then the bijection 0 : Z 2 —y Z 2 given by 


(3.2) 


4 >{x) = x + ei + e 2 , 


which maps Sk to Sk +2 for each k , has the following symmetry property: for all x 
and y , 


(3.3) 


y G Out(x) if and only if <j>(x) G Out(y). 
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FIGURE 9. The doubling graph D, isomorphic to Z, shown on the 
left in correspondence with two successive diagonals 5/, Sk+i of Z 2 . 


Let Dk be the undirected bipartite graph with vertex set Sk U Sk+ 1 , and an edge 
joining x G Sk and y G Sk+i if y E Out (a;). 

The graphs Dk are isomorphic to each other for all fceZ. The doubling graph 
is a generic graph D that is isomorphic to each We can also interpret 1) as 
the image of Z 2 under the equivalence relation x = (j)(x). In the case currently 
under discussion, we can take D to be simply Z, with nearest-neighbour edges, as 
shown in Figure [9j Consider the map v : Z 2 —» Z given by 

(3.4) v((xi,x 2 )) = xi - x 2 . 


Restricted to the set Sk U Sk+ 1 , this gives an isomorphism between Dk and D. for 
any k. 


Recall the definition of the hard-core model as given in Section 1.2 a Gibbs 


distribution for the hard-core model on a graph with vertex set W with activity 
A > 0 is a distribution on configurations rj G {0, 1} W satisfying (1.2). 


Consider the hard-core model on the doubling graph D with vertex set W = Z. 
This is a bipartite graph, with bipartition W = Wo U W\ where Wo and W\ are the 
sets of even and odd integers respectively. We consider the following two update 
procedures for configurations on {0, 1} W . For an “odd” update, for each vertex 
x G W\ independently, resample rj(x) according to the values at its two neighbours, 
setting rj(x) = 0 with probability 1 if either of the neighbours takes value 1, and 
otherwise setting rj(x) = 1 with probability 1 — p. For an “even” update, do the 
same for vertices in Wo- Set A = 1/p—1, so that 1—p = A/(l + A). Since each of 
Wo and W\ is an independent set of D, any Gibbs distribution for the hard-core 
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Figure 10. The Markovian distribution corresponding to a Gibbs 
measure for the hard-core model on the doubling graph W yields 
a Markovian distribution /jl p on each of the two vertex classes Wo 
and W[. Since the Gibbs distribution is invariant under the update 
procedures, the distribution /i v is invariant for the PCA. 


model with activity A is invariant under both of these update operations. (This is 
a version of Glauber dynamics). 

Take some even fceZ. Suppose we start from a configuration on {0, 1} W , which, 
via the isomorphism (3.4) between D and Dp : . under which Wo maps to Sk and W± 
to Sk+i, corresponds to a configuration in {0, l} 5 4 u ‘ s ’/=+i _ Perform an odd update, 
resampling the sites of W\, leading to a new configuration on {0, 1} W . Considering 
now (3.4) as an isomorphism between Dk and Dk- 1 , which maps Wo to Sk and W\ 
to Sk- 1 , the updated configuration on {0,1 } W corresponds to a configuration in 
{0, i} s fc-i u ' s 4, whose values at the sites in Sk are left unchanged. We can interpret 
the update as generating a configuration on Sk -1 from a configuration on Sk- This 
procedure is identical to that which occurs in one iteration of the PCA A p . 

If we then perform an even update, resampling the sites of Wo, we can pass in the 
same way to a configuration on the sites of S’/ c _ 2 U Sk- 1 , which corresponds to the 
next step of the PCA. 


Continuing to perform odd and even updates alternately, we reproduce the evo¬ 
lution of the PCA. A Gibbs distribution on D is characterized by its marginal 
on the vertices of one half of the bipartition, say Wo- Since the distribution is 
preserved by the updates, this distribution on {0, l}^ 0 is 2-periodic for the PCA. 


In fact, for any A there is a unique Gibbs distribution for the hard-core model on 
Z. Since the hard-core interaction is homogeneous and nearest-neighbour, this 
Gibbs distribution is itself a stationary Markov chain indexed by Z. Let Q = Q p 
be its transition matrix. Therefore, the marginal distributions on Wo and W\ are 
in fact equal to each other. Call this marginal distribution /i p . Then /i p is the 
law of the stationary Markov chain with transition matrix P = Q 2 . This /i p is a 
stationary distribution for the PCA A p , and the matrix P is the one in (3.1). See 
Figure 10 for an illustration. 


As mentioned earlier, /i p is a time-reversible stationary distribution for the evolution 
of the PCA. This follows from the reversibility of the process of configurations 
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on G under the update procedure (this is essentially the standard reversibility 
property for Glauber dynamics and is easy to verify by checking the detailed 
balance equations). In fact, from the uniqueness of the Gibbs distribution for the 
hard-core model on G, one can deduce quite easily that there is only one reversible 
stationary distribution for the hard-core PCA A p . However, this argument does 
not preclude the existence of other non-reversible stationary distributions. From 
Theorem [lj we know that such distributions do not in fact exist, but proving this 


required a different argument, as given in Section 2.4 


In contrast, in the next section we will use the implication in the other direction; in 
situations where there exist multiple Gibbs distributions for the hard-core model, 
we can conclude that there are multiple periodic distributions for the corresponding 
PCA; then the PCA is non-ergodic, and draws occur with positive probability in 
the corresponding game. 


3.2. General framework. Recall that in the setting of Theorem [2j we have a 
locally finite graph G with vertex set V, along with a partition (S k : k G Z) of V 


and an integer m > 2 , such that conditions (Al) and (A2) given in Section 1.2 
hold. 


We also defined Dk be the graph with vertex set S k U • • • U Sk+m- 1 , with an 
undirected edge (x, y ) whenever (x, y ) is a (directed) edge of V. For convenience 
we will also use Dk to denote the vertex set Sk U ■ ■ ■ U Sk+ m - i- 


Lemma 3.1. The graphs Dk are isomorphic to each other for all k e Z. 


Proof. Consider the map Xk defined on Dk under which 

Jx if x e S k+ 1 U • • • U Sk+m-1 
|</>(;r) if x E Sk 

From assumptions (Al) and (A2) above, Xk is a graph isomorphism from Dk to 
D^ |-i• Hence indeed Dk and D k + 1 are isomorphic, and so by induction any two 
Dk, D^ are isomorphic. □ 

We then take D to be a graph isomorphic to any Dk . (When m — 2 we sometimes 
call D the doubling graph). Note that D is m-partite. Specifically, let us fix 
some isomorphism / 0 from D 0 to D, and let W t be the image of S) under / 0 , for 
i = 0, ..., m — 1. Then (W 0 ,..., W m _ i ) is a partition of the vertices of D into m 
classes, and assumption (Al) guarantees that there are no edges within a class H 7 *. 

It will be important that we can map both Dk and Dk+i to D in such a way that 
the vertices common to Dk and Dk+\ have the same image in both maps. 
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Lemma 3.2. There exists a family of maps (f k : k G Z) snc/i that fk is a graph 
isomorphism from Dk to D, and such that the following properties hold. 

(a) For each k, 

f ( \ _ I fk+i( x ) f or x E D k n D k+ 1 = S k +i U • • • U S k +m- 1 
\fk+ 1 ( 0 (^)) forxeSk. 

(b) For each k and each r G {fc, k + 1, ..., k + m — 1} 7 the image of S r under 
fk is Wrmodm. 

(c) Let x G S k and y G D k = S k U • • • U S k +m- 1- Then y G Out(x) if and only 
if fk(y) is a neighbour of f k (x) in D. 


Proof. Let / 0 be the isomorphism from D 0 to D described just above. Then we 
can compose / with the isomorphisms \k defined in the proof of Lemma 3.1 , by 
setting 

fo ° X-i ° X -2 ° • • • ° Xk for k < 0 

for k > 0. 


fk = 


fo ° Xo 1 ° Xi 1 ° 


°X fc _i 


Then using assumption (A2), it is easy to check by induction upwards and 
downwards from 0 that f k is an isomorphism from D k to D satisfying the properties 
stated in (a) and (b), for each k. 

Finally note that by (Al), if x G S k then Out(x) C S k U • • • U S k +m- i while In(x) 
is disjoint from Sk U • • • U S k + m - 1 - By definition, the set of neighbours of x in the 
graph D k is then Out(x). Then part (c) follows since f k is a graph isomorphism 
from D k to D. □ 


Given a hard-core configuration in {0,1} D , we can consider Glauber update steps 
that resample the vertices of one of the vertex classes ITo, Wi ,..., W m -\. To 
perform an update of the class Wp for each v G W t independently, let the new 
value at v be 0 if any neighbour of x has value 1, and otherwise let the new value 
at v be 0 with probability p = 1/(1 + A) and 1 with probability 1— p = \/(l + \). 
If a distribution on {0,1} D is a Gibbs distribution for the hard-core model on D 
with activity A = 1/p — 1, then it is invariant under this update procedure for 
each z = 0,1,..., m — 1. (Again, this is a version of the Glauber dynamics for the 
hard-core model on D.) 


Proof of Theorem [j| We start by defining an analogue of the hard-core PCA A p 
in the general setting. As before, we have the recursion (2.1) for the outcome of 
the game started from x G V, in terms of the outcomes started from the elements 
of Out(x) together with the information about whether x itself is open or closed. 
(Recall that we treat the game from x as a win if x is closed.) 
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As in previous sections we can specialize that recursion to configurations involving 
only the symbols 0 — W and 1 = L. This gives the following recursion for a family 
of variables ('y(x) : x G V) G { 0 , 1 } V (which we do not assume to be necessarily 
game outcomes): 

x closed =>• ^{x) = 0; 


(3.5) 


If x G 

(7(x) 


x open 


l{ x ) = 


1 if 7 ( 2 /) = 0 for all y G Out (a;) 

0 if 7 ( 2 /) = 1 for some y G Out(x). 


Ski then Out(x) C Sk+i U • 
x G Sk) in terms of (7(2:) 


U Sk+m-i- Thus, the recursion ( 3 . 5 ) gives 
G Sk- 1-1 U 


x G bk +1 U • • • U Sk+m- 1 ) and the random 
configuration of closed and open sites in Sk (which we take as usual to be product 
measure with each site closed with probability p). This is analogous to the PCA 
A p considered earlier (although for m > 2, a “state” of the PCA is now more 
complicated to describe). 

Fix K G Z, and take some boundary condition (7(0;) : x G Sk U ■ ■ ■ U Sx+m-i): 
which we allow to be be random, but which is independent of the configuration of 
open and closed sites in {J r<K S r . Applying (3.5) repeatedly then generates an 
evolution (7(0;) : x G S r ) r <K+m- i- 

We will couple this evolution with a process of configurations in { 0 , 1 } D . For 
k < K and v G D, define 07(1?) = 7 (/ fc ~ 1 (u))- Then 07 G { 0 , 1 }-° for each k. The 
idea is now to show that the transformation from 07+1 to 07 is identical to a 
hard-core update of the vertex class Pf4 m odm, with randomness provided by the 
configuration of open and closed vertices in Sk- Notice that 07+1 is a function of 
(7(0;) : x G Sk +iU - • -U Sk+m), while cr k is a function of (7(2;) : x £ Sk U - • •US'fc +m _i). 
If v G Wi where i 7^ k mod m, then by Lemma 
Sk+ 1 U • • • U S k+ m- 1- Thus 07+1 (u) = cr k {v). 


3.2 


a) and (b), f k+ 1 (v) = f k (v) G 
the only sites in D which can 
change their value between the configuration 07+1 and the configuration 07 are 
those in Wk modm ■ Consider such a v G Wk mo d mi and let x = /F 1 ^) so that x G Sk 


(by Lemma 3 . 2 (b)) and 07(1?) = 'y(x) 


Translating ( 3 . 5 ) and using Lemma 3 . 2 [ c) gives that for v G Wk mo dm'- 
x closed =>■ ak(v) — 0) 


( 3 . 6 ) 


x open =»■ crfe(u) = 


1 if o’fc+i (u) = 0 for all u ~ v 
0 if 07+1 (w) = 1 for some u ~ v. 


Each bit of randomness (the information about whether x G Sk is open or closed) 
is used only once. Since each x is closed with probability p independently, we have 
that the conditional distribution of ak given 07+1, 07+2,..., ctk is precisely that 
obtained by performing a hard-core update of the vertex class Wkmodm- 
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Now let fi be a Gibbs distribution for the hard-core model on D. By choosing 
the distribution of the boundary condition (7(2;) : x G Sk U • • • U Sx+m-i) 
correspondingly, we can arrange that Ok has distribution /i. But then since ft is 
invariant under the hard-core updates, or has distribution fi for all k < K also. 

So, suppose that there exist multiple Gibbs distributions, and let /j and v be two of 
them. By alternating between ji and u, we can arrange a sequence indexed by K of 
boundary conditions (7^(2;) : x G Sk U ■ • • U Sx+m- 1) that induces a sequence of 
distributions of a q A ' having both p and v as limit points as K —y 00. In particular, 
the configuration does not converge almost surely as K —y 00 (in the product 
topology). So the sequence of configurations (j^(x ) : x G So U • • • U S m _ 1) does 
not converge almost surely. 


Now we apply the same argument that we used for the second part of the proof 
of Proposition T2 If the game started from x is not a draw, then one player 
has a strategy which guarantees a win in fewer than N moves, where N G N 
is an almost surely finite random variable which depends on the percolation 
configuration. Then the values ^ K \x) must agree for all large enough K. Hence 
if there is zero probability of a draw from each site, then the configuration 
(7^(2;) : x G So U • • ■ U S m - 1) converges almost surely as K —> 00. By the 
argument in the previous paragraph, this contradicts the existence of multiple 
hard-core Gibbs distributions. □ 


3.3. Remarks on the converse direction. We do not know whether the con¬ 
verse of Theorem [2] holds; that is, whether the uniqueness of the hard-core Gibbs 
distribution on D implies that there are no draws on G. Two separate issues arise. 


First, assume that the hard-core Gibbs distribution is unique, and suppose we 
wish to conclude that th ere is a unique distribution of the evolution (7(2;) : x G 
G) G { 0 , 1 } V satisfying ( 3 . 5 ) that is stationary in the sense that (7(2;) :iGG) 
and (7(0(2;) :igG) have the same law. Uniqueness of the hard-core distribution 
gives only that there is a unique such law satisfying an additional condition, which 
is that ay has the same distribution for all k (and not merely for any two k 
which agree modm). In the case m = 2, this condition can be interpreted as a 
reversibility property, although the precise meaning of this is complicated by the 
lack of a canonical bijection between Sk and Sk+i- Specifically, the condition is 
equivalent to the statement that the ordered pair [(7(2;) : x G So), (7(2/) : V 7 S_i)] 
has the same distribution as [(7(2;) : x G S 0 ), (7 4 >{ y ) '■ y € S_i)]. In the example 
of the game in Z 2 , this is also equivalent to symmetry of the law of the evolution 
on Z 2 under reflection in the line 27 + X2 = 0 . 
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However, there may be other stationary evolutions that do not satisfy such a 
reversibility condition. (As mentioned before, our argument for Z 2 rules these out 
by a different method.) 


Second, note that when G is bipartite (for example, whenever m — 2), the 
uniqueness of law of a stationary evolution with the binary alphabet {0,1} is 
enough to imply that the probability of draws is zero. This follows by an argument 
similar to the proof of Proposition 2.1 showing that the PCA on three symbols 
{0, ?, 1} is ergodic whenever the binary PCA is ergodic. 

However, if G is not bipartite (so that there are sites which can be reached on 
the turn of either one of the players), then this argument breaks down. The 
monotonicity associated to the stochastic order generated by 0 < ? > 1 (see 
Lemma 2.1 (ii)) still applies, but the monotonicity associated to the stochastic 
order generated by 0 < ? < 1 (see Lemma 2.1 4)) fails; the operation analogous to 
a single application of the operator F p involves moving from a configuration on 
the vertices Sk+i U • • • U Sk+ m to a configuration on the vertices Sk U • • • U Sk+ m - i, 
and these two sets overlap when m > 2. 


Hence it is no longer clear that uniqueness for the binary PCA would imply the 
same for the three-symbol PCA, as required to conclude that there are no draws. 
Closely related examples in which ergodicity of a binary PCA does not imply 
ergodicity of its envelope PCA are noted in [ BMM13] , 


3.4. Example graphs with d > 3. We now give several examples of graphs G 
to which one may hope to apply Theorem |2j (As we will see, the result can indeed 
be applied in some cases, but not in others.) We consider the vertex set V = Z d 
(or subsets thereof), for various different choices of the set Out (a:) of vertices to 
which the token can move from x. In each case we consider the percolation game 
in which each site is closed with probability p and open with probability 1 — p. 
All our examples can be regarded as natural extensions of the original Z 2 game, 
in the sense that they reduce to it when we set d = 2. 

Example 3.1. Let Out (a?) = {x + e* : 1 < i < d}. So | Out(x)| = d. This is 
perhaps the most natural extension of all. However we cannot apply Theorem 
[2] because there is no choice of the automorphism </> for which assumption (A2) 
holds. 


Example 3.2. Let Out(x) = {x ± : 1 < i < d — 1}. This is the example 

already mentioned in Section 1.2 Here |Out(x)| = 2 (d — 1). 


preserves parity, it is natural to restrict to the set of even sites Z; 
Y^, x i is even}. 


Since any step 

d _ f ^ 'T/d 

even 


:= {x e 
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FIGURE 11. The three-dimensional body-centred cubic lattice 
(which is the doubling graph for the PCA associated to Exam¬ 
ple 3.3 when d — 4). The two underlying copies of Z 3 are shown 
with red and with blue vertices. The black lines are the edges of 
the body-centred cubic lattice, and the dotted red lines show the 
nearest-neighbour edges in the red copy of Z 3 . 


In two dimensions, the game is isomorphic to the original game on Z 2 . For general 
d, conditions (Al) and (A 2 ) hold with m = 2 if we set S k = {i£ Zf ven '■ x d = k} 
and 4>{x) = x + 2 e d . 

To obtain the doubling graph, consider D k = S k U S k+ 1 with an edge between 
x G S k and y G S k + 1 whenever y G Out(x). This gives a graph isomorphic to the 
standard cubic lattice Z ^” 1 (for example, the map 

(a?i, — , x d _ i, x d ) (xi,..., x d _ i) 

gives a graph isomorphism from D 0 = S 0 U Si to Z d_1 ). 

The graph is vertex-transitive. By Theorem [2j if there exist multiple Gibbs 
distributions for the hard-core model on Z d_1 with activity A, then the percolation 
game on G with p = 1/(1 + A) has positive probability of a draw from any vertex. 


Example 3.3. Now let Out(a:) = {x±ei ±e 2 • • • ±e,j-i + e d }, so that | Out(x)| = 
2 d ~ 1 . Each step changes the parity of every coordinate, so we restrict to the 
set Z^ cc = {x G Z d : Xi = x 3 mod 2 for all i,j}- Putting an edge between x 
and y whenever y G Out(x), we obtain the body-centred cubic lattice in d 
dimensions. This consists of two copies of (2Z) d , each offset from the other by 
( 1 , 1 ,..., 1 ), so that each point of one lies at the centre of a unit cube of the other; 
the edges are given by joining each point to the 2 d corners of the surrounding unit 
cube. See Figure 11 for an illustration. 

Conditions (Al) and (A2) hold for m = 2 with S k = {x G Z^ cc : x d = k} 
and cj)(x) — x + 2e d . The doubling graph D isomorphic to D k = S k U S k +i for 
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Figure 12. Triangular lattice and hexagonal lattice. 


each k is now the body-centered cubic lattice in d — 1 dimensions. The map 
v(x) = (xi,x 2 , ■ • • ,Xd~ i) from Z^ cc to restricts to an isomorphism between 
Dj. and D for each k. 


When d = 2 or d = 3 the graph G is isomorphic to that in Example 3.2 above, but 
for d > 4 the graphs are different. Existence of multiple hard-core distributions 
on will imply existence of draws on Z^ cc . 


Example 3.4. Let Out(x) = {# + 22ies e* : S C {1, ..., d} with 1 < \S\ < d — 1 }. 
So a move of the game corresponds to incrementing at least one, but not all, of 
the coordinates by one. Then | Out (a?) = 2 d — 2. 

Conditions (Al) and (A2) hold with m = d if we set Sk — {x : 22 x i — and 
4>(x) = x + e\ + e 2 + • • • + ed- For d — 2 the game is the same as ever. For d > 2 
there are some new features. For the first time we have m > 2, and the graph G is 
not bipartite; from a given starting vertex, there are vertices that can be reached 
when it is either player’s turn. The graph D is (d — l)-dimensional and d-partite. 
For d = 3, it corresponds to the triangular lattice. For example, the map 

(3.7) (xi,x 2 ,x 3 ) —>■ xi(l,0) +x 2 {-\, +£ 3 (-7,-^p) 

is an isomorphism from Dk = Sk U • • • U Sk+m-i to the triangular lattice for each 
k. 


Example 3.5. Fix r with 1 < r < d, and now restrict to sites x G such 
that 22 x i = 0 or r mod d. For x with 22 x i = 0 mod d, let Out (a;) = {x + 
22 i&s O, f° r any S C {1,... ,d} with IS) = r}. Meanwhile, for x with 22 x i = 
r mod d, let Out(x) = {a: + 22i^s e u f° r any S C {1,..., d} with \S\ = d — r}. 

Now | Out (a:) | = (^) for all x. Replacing r by d — r gives an isomorphic graph, so 
we may assume 1 < r < d/2. Then conditions (Al) and (A2) hold with m = 2, 
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FIGURE 13. The cubic lattice and the diamond cubic graph. 


with (fi(x) = x + Y^%=\ e ii an( i with Sk = {x : ^2 Xi = dk/ 2 } for even k and 
Sk = {x : Yh x i — d{k — l )/2 + r} for odd k. 

For d — 2 (and hence r = 1) the game is the familiar two-dimensional game. For 
d — 3 and r = 1, we get | Out(x)| = 3 and the doubling graph D is the two- 
dimensional hexagonal lattice; this is the image of {x G Z d : Y2 x i = 0 or 1 mod 3}, 
with edges between x and y where y G Out(x), under the map (3.7) above. 

For d — 4 and r = 2, the graph G is isomorphic to the d — 4 case of Example |3.2 
above, and so D is the standard cubic lattice Z 3 . For d — 4 and r = 1, we have 
| Out (a:) | = 4, and D is the so-called diamond cubic graph (see for example 
Section 6.4 of |CS99j ). This graph may, for example, be represented as 

{y G Z 3 : y 1 = y 2 = y 3 mod 2 and y x + y 2 + y 3 = 0 or 1 mod 4}, 

with edges between nearest neighbours (which are at distance a/ 3/4). This is the 
image of {x G Z d : Y2 x i = 0 or 1 mod 3}, with edges between x and y where 
y G Out(x), under the map 

x i — x 2 - 

-Xi + X 2 - 
Xi + X 2 - 

(see | NS08] ). See Figure 13 for an illustration. 




Theorem 3. There is positive probability of a draw from every vertex for suffi¬ 
ciently small p in the following cases: Example T2 for all d > 3; Example \3.3\ for 
d = 3 and d = 4; Example 3.4 for d = 3; and Example 3.5 for d = 3 (with r = 1) 


and d = 4 (with r = 1 or r = 2). 
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Proof. In the cases listed, it is known that there exist multiple Gibbs distributions 
for the hard-core model on the associated graph D when the activity parameter A 
is sufficiently high. For the standard cubic lattice in any dimension greater than 
1, the result goes back to Dobrushin [Dobfiol . Other models in two and three 
dimensions were covered by Heilmann |Hei74] and Runnels [Run 75] , including the 
triangular and hexagonal lattices in two dimensions and the body-centered cubic 
lattice and the diamond cubic graph in three dimensions. 

Theorem [2] shows that there is positive probability of a draw for small enough p 
from some vertex, and since all the graphs G are vertex-transitive, the conclusion 
holds for every vertex. □ 


It is expected that in fact the hard-core model on D has multiple Gibbs distributions 
for A sufficiently large in all of Examples 3.2 3.5 whenever d > 3 (so that D has 
dimension at least 2). This could likely be proved by Peierls contour arguments, 
although this requires a suitable definition of a contour, which is typically graph- 
dependent, and less straightforward than in other settings such as the Ising model. 
Via Theorem [2j such non-uniqueness would imply existence of draws for the 
corresponding graphs G. 

We emphasize again that there is a more fundamental ob stacl e to proving existence 
of draws for the standard oriented lattice of Example 3T in that our condition 
(A2) does not hold here. 


3.5. Extending the hard-core correspondence. Various further extensions 
can be made while still preserving the correspondence to the hard-core model. For 
example, in the class of models considered in Theorem [2j we can augment the set 
of allowable moves from site x to include the point <f>{x) itself. 

Specifically, replace (Al) and (A2) by the following assumptions: 

(AT) For all x G Sk, Out(x) C S^+i U • • • U Sk+ m - 1 U 

(A2 7 ) There is a graph automorphism <f of G that maps Sk to Sk+m f° r every 
k, such that (p(x) G Out(x) for all x, and such that Out (a;) \ {0(x)} = 
ln(0(x)) \ x. 

Define D as before; D is a graph isomorphic to any D^, where D is the graph 
with vertex set Sk U ■ ■ • U Sk+ m - 1 and an undirected edge (x, y ) whenever (x, y ) is a 
directed edge of V. Note now an edge (x, </>(x)) in G does not give a corresponding 
edge in D. 

Proposition 3.1. Suppose that the graph G satisfies (AV ) and (A2'). If there 
are multiple Gibbs distributions for the hard-core model on D with activity A < 1. 
then the percolation game on G with p = 1 — A has positive probability of a draw 
from some vertex. 












32 


HOLROYD, MARCOVICI, AND MARTIN 


The method of proof is a slight variation on that of Theorem [2| which we indicate 
briefly. To reflect the presence of the edge (x,(f>(x)), we change the hard-core 
update procedure. When we perform an update of the vertex class W t . we now 
add that any vertex v G W t which is in state 1 before the update must move to 
state 0 after the update; otherwise the update at v proceeds as before. 

Again one can show that hard-core Gibbs distributions are stationarity under 
such updates, but now with the activity parameter A equal to 1 — p rather than 
1/p — 1 as previously. (To verify the stationarity, one can start by checking the 
detailed balance condition for an update at a single site; then if the distribution is 
stationary for the update at any single site, it is also invariant under simultaneous 
updates at any set of non-neighbouring sites.) 


Note that now when p —y 0, we have A —> 1 rather than A —>• oo. Hence to show 
existence of draws for some p, we need multiplicity of Gibbs distributions for some 
A < 1. For the case of the standard cubic lattice, Galvin and Kahn |GK04 ] show 
that this holds for sufficiently high dimension, so that we can deduce the existence 
of draws for the variant of Example 3.2 in which Out(x) = {x ± e* + : 1 < % < 

d — 1} U {x + e^}, when d is sufficiently large. 


4. Open questions 

4.1. The oriented cubic lattice. For the percolation game on Z d (where the 
allowed moves from site x are to any open site x + Ci for 1 < i < d ), do there exist 
any p G (0,1) and d > 3 for which draws occur with positive probability? 

4.2. Monotonicity and phase transition. For lattices where draws are known 

to occur (for example, the even sublattice of 7L d with d > 3 and moves allowed 
from x to any open x + ± e* for 1 < i < d — 1), is the probability of a draw 

starting from the origin non-increasing in the density p of closed sites? Or, at 
least, is the set of p that have positive draw probability a single interval containing 
0 (so that there is a single critical point at the upper end of the interval)? If so, 
what happens at the critical point? (It is also unknown whether the hard core 
model on has a single critical point for uniqueness of Gibbs distributions; see 
e.g. [ BGRT13 ] for discussion and recent bounds.) 


4.3. Reversibility. For lattices where draws are known to occur, such as the 
example described in question 4.2 above, do there exist {0, l}-valued invariant 
distributions for the game evolution that are not projections of a hard-core Gibbs 
distribution on the graph I), and therefore lack the corresponding reversibility 
property? Or is it the case that the game evolution has multiple invariant 
distributions precisely when there are multiple hard-core Gibbs distributions? 
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4.4. Target game win probability. Compute the winning probability for the 
target game on Z 2 (in which the winner is the first player to move to a closed 
site). 

4.5. Misere games. Besides our target game, there are other natural misere 
variants. For example, suppose as before that moves are allowed from x G Z 2 
to an open site in { x + e\,x + e? }, but now declare that if both these sites are 
closed, then the player whose turn it is to move wins. Does this game have zero 
probability of a draw for all pi Unlike the original percolation game or the target 
game, there is apparently no useful correspondence to a PCA with alphabet {0,1}. 

4.6. Elementary probabilistic cellular automata. Is every elementary PCA 
(i.e. one with 2 states and a size-2 neighhborhood) on Z with positive rates ergodic? 
Can our weighting approach be extended to prove ergodicity for other PCA in 
this class? 

4.7. Undirected lattices. The following game is considered in [ BHMW15 ]. As 
usual, each site of 7L d is independently closed with probability p, and two players 
alternately move a token. From an open site x, a move is permitted to any open 
nearest neighbour x ± e* provided it has not been visited previously. A player who 
cannot move loses. This game is closely related to maximum matchings, and this 
is used in [ BHMW15 ] to derive results for biased variants in which odd and even 
sites have different percolation parameters. However, for the unbiased version 
described above it is unknown whether there exist any p > 0 and d > 2 for which 
draws occur with positive probability. 
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